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Target  Detection  Using  a  Linear  Sum  of  Matched  Filter  Outputs 


1.  INTRODUCTION 

The  matched  filter  is  the  most  common  receiver  structure  used  for  target  detection  in  echolo- 
cation  systems.  As  shown  in  Figure  1,  it  is  composed  of  a  multiplier  that  combines  the  received 
echo,  5(t),  and  the  processing  waveform,  g{t  —  r^),  where  Td  is  the  range  delay.  This  is  followed  by 
an  integrator.  Since  the  processing  signal  is  a  function  of  range  delay ^ ,  the  complex  matched  filter 
output,  y,  is  also  a  function  of  range  delay  and  is  sampled  in  time  to  form  “bins.”  The  amplitude 
of  the  resulting  complex  value  in  each  bin,  the  “detection  statistic,”  is  compared  to  a  threshold. 
When  the  statistic  exceeds  the  threshold,  a  detection  is  declared. 

The  above  detection  scheme  is  based  on  only  one  bin  value  at  a  fixed  range  being  available 
during  a  single  ping  cycle,  and  is  thus  referred  to  as  “single-observation  detection.”  In  some 
cases,  multiple  observations  are  available  because,  for  example,  an  echo  may  be  divided  into  sub¬ 
bands  (in  frequency)  that  are  separately  processed  (matched  filtered).  When  several  observations 
are  simultaneously  available,  they  can  (or  should)  be  used  collectively  to  determine  if  a  target  is 
present.  This  can  be  accomplished  by  using  a  more  general  form  of  the  matched  filter  where  the 
observations,  each  originating  from  a  separate  matched  filter,  are  summed  to  form  a  single  detection 
statistic. 

This  report  documents  mathematical  techniques  for  approximating  the  probability  density 
function  (PDF)  and  other  statistical  functions  that  characterize  matched  filter  outputs  and  detec¬ 
tion  statistics  that  are  a  linear  combination  of  the  outputs.  The  approximations  can  be  used  for 
computing  detection  thresholds  and  detection  probabilities  a  function  of  signal-to-clutter  ratio 
(SCR).  The  approximations  are  either  based  on  Fade  approximations  to  the  multi-dimensional 
characteristic  function  (CF),  or  are  discrete  mixtures  based  on  an  exact  or  approximate  form  of 
the  Rayleigh  mixture  parameter  (echo  power)  PDF. 

This  report  documents  mathematical  techniques  for  determining  the  probability  of  detection, 
Prf,  and  detection  threshold,  7,  for  a  given  probability  of  false  alarm,  when  using  either  a 
single  output  or  a  linear  sum  of  multiple  outputs  of  a  matched  filter.  This  includes  techniques  for 
modeling  the  statistical  nature  of  the  matched  filter  response  to  target  and  clutter  echoes.  Specific 
results  include  the  following: 

•  Recently  reported  techniques  for  approximating  the  statistical  functions  that  model  a  single 
matched  filter  output  and  are  based  on  Fade  approximations  to  the  CF  are  extended  to  multiple 
observations  (Section  3.5). 

•  Techniques  were  developed  that  yield  accurate,  discrete  Rayleigh  mixture  models  of  multidi¬ 
mensional  PDFs  and  related  statistical  functions  (Sections  2.5  and  3.6). 

Manuscript  approved  June  16,  1999, 

*More  generally,  the  processing  signal  can  also  be  a  function  of  other  variables  such  as  Doppler  and  bearing. 
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8*(t-T^) 

Fig.  1  -  The  matched  filter. 


•  The  aforementioned  Pade-based  techniques  have  been  extended  to  include  approximations  to 
the  statistical  functions  that  model  the  Rayleigh  mixture  parameter  (Section  2.4). 

•  It  was  found  that  Pade-based  approximations  are  best  suited  for  approximating  the  “tails”  of 
the  PDFs,  and  thus  may  be  used  to  determine  the  detection  threshold  given  a  specific  false 
alarm  probability  where  typically  P^a  <0.001  (Section  3.5). 

•  It  was  found  that  approximations  based  on  mixture  models  are  best  suited  for  approximating 

the  “peaks”  of  the  PDFs,  and  thus  may  be  used  to  determine  target  detection  probabilities 
where  typically  0.1  <  <  0.9  (Sections  2.5  and  3.10). 

Unlike  other  treatments  of  the  detection  problem  where  the  target  is  assumed  to  be  “steady,” 
it  is  assumed  here  that  a  matched  filter  output  fluctuates  in  response  to  the  target  echo.  This  is 
true  in  echolocation  systems  (especially  sonar  systems  operating  in  “search  mode”)  where  the  range 
resolution  is  equal  to  or  larger  than  the  line-of-sight  target  size. 

To  summarize,  this  report  is  organized  into  two  main  sections:  (i)  modeling  a  single  observa¬ 
tion  and  (ii)  modeling  multiple  observations.  The  first  section  provides  the  necessary  background 
from  which  detection  using  multiple  observations  is  developed.  It  begins  with  a  summary  of  the 
underlying  assumptions  of  the  single-observation  detection  problem  for  echolocation  systems.  This 
leads  to  a  discussion  of  the  CF  and  Rayleigh  mixture  model  representations  and  their  correspond¬ 
ing  approximations  for  the  matched  filter  output  statistical  functions.  Examples  of  both  Rayleigh 
and  non-Rayleigh  fluctuation  models  are  used  to  illustrate  the  techniques.  The  second  section  then 
builds  on  these  results  to  account  for  multiple  observations  by  developing  techniques  to  approxi¬ 
mate  multivariate  statistical  functions.  It  includes  a  description  of  the  receiver  structure  and  the 
associated  calculations  of  the  detection  threshold  and  probability  of  detection.  Examples  of  receiver 
operating  characteristic  curves  conclude  the  section. 

2.  BACKGROUND:  MODELING  A  SINGLE  OBSERVATION 

2.1  Statistical  Assumptions 

Suppose  the  complex  output  of  a  matched  filter  receiver  is  given  by 

y  zr:  u-\-  jv  =  (1) 

where  r  and  9  are  respectively  the  amplitude  and  phase  of  the  matched  filter  output.  Under  the 
common  assumption  that  there  is  no  knowledge  of  its  likely  value,  the  phase  is  modeled  as  a  uniform 
random  variable  in  the  interval  [0,27r).  This  implies  that  u  and  v  are  zero-mean  random  variables 
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and  their  joint  PDF  is  circularly  symmetric.  Using  the  polar  notation  (r,^),  this  joint  PDF  has 
the  following  form: 


where  /^(r)  is  the  PDF  of  the  amplitude  r.  The  CF  for  u  and  v  is  equal  to  the  two-dimensional 
Fourier  transform  of  the  joint  PDF  in  (2)  and  is  denoted  as  It  is  possible  to  show  that 

due  to  the  circular  symmetry  of  the  joint  PDF,  the  CF  is  also  circularly  symmetric  and  is  only  a 
function  of  the  “radial  frequency,”  p  =  -f-  Consequently,  only  the  relationship  between  the 

CF  and  the  amplitude  PDF  is  important  and  is  expressed  by  the  following  Hankel  transform  pair 

[10]: 

roo 

Hp)=  /  fR{'>')Jo{pr)dr  (3) 

Jo 

roo 

fnir)  =  r  p^(p)Jo(pr)  dp  (4) 

Jo 

where  Jo(*)  is  the  zeroth-order  Bessel  function  of  the  first  kind.  Detection  probabilities  are  given 
by  the  survival  function  (SF),  defined  as 

/oo 

fR(t)dt.  (5) 


2.2  The  Rayleigh  Mixture  Amplitude  Fluctuation  Model 
The  Rayleigh  mixture  amplitude  fluctuation  model  is  defined  as 

7  (“  y) 

where  /t('7‘)  is  the  PDF  of  the  Rayleigh  parameter,  r,  and  is  referred  to  as  the  "characteristic 
PDF.”  Integration  and  Hankel  transformation  of  (6)  leads  to 


Snir)  =  jT  exp  ( 

(7) 

$(/>)  =  exp  I 

yi)MrUr. 

(8) 

This  model  was  deduced  from  the  observation  that  some  radar  clutter  appears  to  exhibit  Rayleigh 
statistics  whose  power  is  dependent  upon  the  portion  of  the  environment  being  observed  [16].  One 
PDF  that  has  been  successfully  employed  as  a  radar  clutter  model  is  the  K-type  PDF,  which  can 
be  shown  to  be  a  Rayleigh  mixture  whose  characteristic  PDF  is  a  gamma  distribution.  Anecdotal 
evidence  suggests  that  acoustic  scattering  from  the  ocean  floor  may  also  follow  a  K-type  PDF,  thus 
suggesting  that  Rayleigh  mixture  models  might  be  employed  as  sonar  reverberation  models^. 


^This  is  a  result  of  data  analysis  performed  by  the  authors.  References  on  this  work  are  not  currently  available. 
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2.3  Approximate  Amplitude  Fluctuation  Models  Based  on  Fade  Approximations 


It  is  shown  in  References  5  and  6  that  the  CF  may  be  approximated  by  a  Fade  approximation 
(rational  approximation); 


$(p)  ss 


1  +  Uip^  +  a2P*  +  ■  •  • 


(9) 


where  L  <  M.  The  coefficients  of  the  numerator  and  denominator  of  (9)  can  be  found  from  the 
amplitude  PDF  by  either  “moment  matching”  at  the  origin  (p  =  0)  or  by  “interpolation.”  In 
the  case  of  moment  matching,  the  denominator  coefficients  are  found  first  by  solving  the  matrix 
equation 


^  ^L-M+l  Cl-M+2 
^L-M+2  f^L-M+3 

where  Cq  =  Ij  Cfc  =  0  for  A:  <  0  and 


\ 

bu-i 

— 

Cz,+M-1  / 

br  J 

Ck 


(A:!)2  22'= 


for  A;  >  0 


(10) 


(11) 


with  denoting  the  2A:-th  moment  of  the  amplitude  PDF.  Once  6i, . .  .,6m  have  been  found, 

the  numerator  coefficients  can  be  determined  using 


«-i 

«<  =  6i +Ci +  ^6i_fcCjt  for  i  =  1,...,X. 

fc=i 


The  Fade  approximation  in  (9)  can  be  written  as  a  partial  fraction  expansion: 


(12) 


$(p) « 


A, 

(p^-p,)'"*'  ■ 


(13) 


Treating  the  CF  as  a  Hankel  transform  and  inverting  the  partial  fraction  expansion  in  (13)  yields 
an  approximation  for  the  amplitude  PDF: 


Integration  of  (14)  yields  an  approximation  for  the  SF; 

c  2Afcr*"'=/i:„,^(y-p^  r) 


(14) 


(15) 
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2-4  Approximate  Rayleigh  Parameter  Fluctuation  Models  Based  on 
Fade  Approximations 

A  Fade  approximation  to  the  CF  can  also  be  used  to  derive  an  approximation  to  the  charac¬ 
teristic  PDF  and  its  associated  SF.  These  approximations  are 


(16) 

(17) 

where  r(n,  t)  denotes  the  incomplete  gamma  function.  The  approximation  in  (16)  may  be  verified  by 
substituting  it  into  (8)  and  carrying  out  the  elementary  integration  to  get  (13).  The  approximation 
in  (17)  follows  from  the  integration  of  (16).  If  mjt  =  1  for  all  k  (which  is  usually  the  case),  then 
(17) reduces  to 

^  (^)  * 

These  equations  are  valid  only  when  Re{pfc}  <  0  for  all  k,  since  this  condition  implies  that  the  terms 
in  (16)  and  (17)  will  eventually  decay  to  zero  as  the  Rayleigh  parameter  becomes  large^.  In  the 
following  examples,  the  Nakagami,  K-type  and  Rayleigh  fluctuation  models  are  used  to  illustrate 
the  accuracy  of  (16)  and  (17)  by  comparing  them  to  their  known  analytic  forms. 


The  Nakagami  Target  Fluctuation  Model 


Consider  the  example  of  a  Nakagami  target  fluctuation  model  given  by: 


/«(»■)  = 


1  /  Nr'^\ 

,  ^(-i)'=iV(Ar-n)...(Ar  +  fc-i)/uVV 
=  i  +  ^ - (M)5 - yw)  • 


(19) 

<20) 

(21) 


This  fluctuation  model  follows  from  the  chi-square  model  for  target  strength  (r®)  [12,7].  In  par¬ 
ticular,  JV  =  1  corresponds  to  a  Rayleigh  model  and  N  =  2  corresponds  to  a  one-dominant-plus- 
Rayleigh  model  [4].  For  A  <  1,  the  characteristic  PDF  for  this  model  is  given  by  [14] 


frir) 


Ny  ~  (-Nr/a^y 

r(A)^r(i-i-i)r(i-A-i) 


for  0  <  r  <  a^/N 


I  0 


for  a^/N  <  T. 


(22) 


^It  is  shown  in  Reference  5  that  for  the  Hankel  transform  to  be  “stable,’’  which  results  in  approximations  to  the 
amplitude  PDF  and  SF  that  decay  to  zero  as  r  — oo,  the  poles  (p^  for  A:  =  1 . . ,  M)  cannot  be  positive  real  numbers. 
This  admits  the  possibility  of  complex  poles  with  positive  real  parts  and  non-zero  imaginary  parts.  Such  poles  would 
not  yield  approximations  to  the  characteristic  PDF  and  SF  that  decay  to  zero  as  r  — ►  oo.  Nevertheless,  it  has  been 
the  authors’  experience  that  poles  for  Pade-based  approximations  to  the  amplitude  PDFs  and  SFs  for  models  that 
are  known  to  be  Rayleigh  mixtures  (e.g.,  the  Weibull,  Nakagami  and  K-type  models)  have  negative  real  parts. 
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The  exact  form  of  the  characteristic  SF  is  obtained  by  integration  of  (22)  and  is  given  by 

r  /'NY  Y  ~  (-Nr/a^y 

r(iV)  £  (iV  +  Or(i  +  l)r(l  -N-i)  forO<r<a/iV 


^r(r)  =  \ 

^  0  for  a'^/N  <  t. 

For  N  =  0.5  and  cr^  =  1,  it  is  shown  in  Reference  5  that 

,  -19.78778 +  J258.7502  -19.78778  -  j258.7502 

+  7.782720  +  jl.618012  +  7.782720  -  jl.618012 

30.66573  -  ;102.2120  30.66573  +  ^102.2120 

+  6.813219  +  J4.909758  p"  +  6.813219  -  j4.909758 
-8.790265  +  J9.696267  -8.790265  -  j9.696267 

"^p^  +  4.478996  +  j8.421761  p^  +  4.478996  -  j8.421761  ’ 


(23) 


(24) 


Substituting  the  coefficients  and  poles  of  the  partial  fraction  expansion  in  (24)  into  (16)  and  (17) 
and  carrying  out  the  tedious  algebra  yields 

/r(r)  «  259.5057  exp  (-3.891360r)  cos  (0.809006r  -  1.647122) 

+106.7131  exp  (-3.406610r)  cos  (2.454879r  +  1.279320) 

+13.08764  exp  (-2.239498r)  cos  (4.210881r  -  2.307225)  (25) 

Srir)  «  65.29159  exp  (-3.891360r)  cos  (0.809006r  -  1.442144) 

+25.41406  exp  (-3.406610r)  cos  (2.454879r  +  1.903753) 

+2.744104  exp  (-2.239498r)  cos  (4.210881r  -  1.225219) .  (26) 


The  exact  forms  and  approximations  to  the  characteristic  PDF  and  SF  are  shown  in  Figures 
2  and  3.  Figure  3  reveal  that  (26)  is  a  respectable  approximation  to  the  characteristic  SF  that 
deviates  about  the  exact  form.  This  modeling  error  is  reflected  by  the  more  obvious  deviations 
seen  in  the  approximation  to  the  characteristic  PDF.  Figure  4  shows  the  relative  errors  between 
the  exact  form  of  the  characteristic  PDF  and  SF  and  their  approximations,  and  are  defined  by 


RE  =  logio 


1- 


approximate  equation 
exact  equation 


(27) 


The  negative  value  of  the  relative  error  is  approximately  equal  to  the  number  of  significant  digits 
of  agreement  between  the  exact  equation  and  its  approximation.  Generally,  an  approximation  is 
considered  good  when  RE  <  —2.  For  Figure  4  the  relative  error  was  not  calculated  for  r  >  2 
because  the  exact  form  is  identically  zero  for  this  interval. 


The  K-type  Clutter  Fluctuation  Model 


As  another  example,  consider  the  K-type  clutter  fluctuation  model  given  by 


^'W=r(.-n)(2) 

v+1 

1  K,{br) 

*^+1 

1  K^+i  (br) 

(28) 

(29) 

(30) 
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Fig.  2  -  The  exact  and  approximate  forms  of  the  characteristic  PDF  for  the  Nakagami  target  fluctuation  model. 
The  approximation  is  based  on  the  Fade  approximation  to  the  amplitude  CF. 


Rayleigh  Parameter  (t) 


Fig.  3  ~  The  exact  and  approximate  forms  of  the  characteristic  SF  for  the  Nakagami  target  fluctuation  model. 
The  approximation  is  based  on  the  Fade  approximation  to  the  amplitude  CF. 
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Fig.  4  -  The  relative  errors  for  the  approximations  to  the  characteristic  SF  and  PDF  for  the  Nakagami  target 
fluctuation  model  that  are  shown  in  the  previous  figures.  For  r  >  2,  the  relative  error  is  infinite. 


It  is  weU  known  that  the  characteristic  PDF  and  SF  for  this  model  are  given  by  [14] 


^  .  s  r  (i/  + 1 , 6V/2) 

-  “  T(u  +  1)  • 


(31) 

(32) 


For  1/  =  0.5  and  b  —  -y/S,  it  is  shown  in  Reference  5  that 


$(/))« 


10.04852 
+  2.932845 


7.166959  1.600556  0.9000533 

+  3.403857  p^  +  5.813469  p^  +  19.84990  ‘ 


(33) 


Substituting  the  coefficients  and  poles  of  the  partial  fraction  expansion  in  (33)  into  (16)  and  (17) 
and  carrying  out  the  algebra  yields 


/^(r)  Ri  -0.4500269  exp (-9.924915r)  -  0.8002768  exp (-2.906733r) 

-3.583477  exp(-1.7019287r)  +  5.024257  exp(-1.4664224T)  (34) 

Sr{T)  -0.04534315  exp  (-9.924915r)  -  0.27531830  exp  (-2.906733r) 

-2.1055389exp(-1.7019287r)  +  3.4262004exp(-1.4664224r).  (35) 


These  approximations  cannot  be  visually  distinguished  from  their  comparable  exact  forms  in  (31) 
and  (32)  when  displayed  on  the  same  graph.  Their  differences  can  only  be  seen  from  examining 
the  relative  error  curves  shown  in  Figure  5.  They  reveal  that  the  approximations  generally  agree 
with  the  exact  forms  to  at  least  2  decimal  places. 
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Fig.  5  -  The  relative  errors  for  the  approximations  to  the  characteristic  SF  and  PDF  for 
the  K-type  clutter  fluctuation  model. 


The  Rayleigh  Fluctuation  Model 

As  a  final  example,  consider  the  Rayleigh  target  and  clutter  fluctuation  model  given  by 


/RW  =  iexp(-^)  (36) 

5fi(r)  =  exp^-^^  (37) 

$(/!))  =  exp 

For  =  1  it  is  easy  to  see  that 

f^(T)  =  S(T-<7^)  (39) 

Sr(T)  =  u,  -  r)  (40) 


where  «,(•)  is  the  unit  step  function  and  S(-)  is  its  derivative,  the  Dirac  delta  function.  In  Reference 
5  it  is  shown  that 


-204.5260  +  J47.68627  -204.5260  -  j47.68627 

^  p^  +  10.97484  +  ;5.196570  +  10.97484  -  i5.196570 

58.07762- J28.55325  ,  58.07762  + j28.55325 

/)2  + 8.291362 +  J10.38788  +  8.291362  -  jlO.38788 

-4.878310  + J3.723561  ,  -4.878310  -  j3.723561 

p2  ^  2.839304  +  J15.61901  p^  +  2.839304  -  jl5.61901 
302.6534 
p2  ^  11.80429  ■ 


(41) 


Again,  substituting  the  coefficients  and  poles  of  the  partial  fraction  expansion  in  (41)  into  (16)  and 
(17)  and  carrying  out  the  algebra  yields 
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Fig.  6  -  The  exact  and  approximate  forms  of  the  characteristic  PDF  for  the  Rayleigh  clutter  fluctuation  model. 
The  exact  form  is  a  Dirac  delta  function  with  unity  weight.  The  approximation  is  based  on  the  Fade  approximation 
to  the  amplitude  CF. 


/r(r)  PS!  210.0116  exp(-5.487420r)  cos(2.598285r  -  2.912529803) 

+64.71706  exp(-4.145681r)  cos(5.193940r  +  0.456936839) 
+6.137004exp(-1.419652T)  cos(7.809505r  -  2.489640610) 

+151.3267  exp(-5.902145r)  (42) 

Sr{T)  «  34.58986  exp  (-5.487420r)  cos  (2.598285r-  2.470307) 
+9.738367exp(-4.145681r)cos(5.193940r+  1.354105) 

+0.7731666  exp  (-1.419652r)  cos  (7.809505r  -  1.09867) 

+25.63927  exp  (-5.902145r).  (43) 

Figures  6  and  7  reveal  that  (42)  and  (43)  are  crude  approximations  to  a  delta  function  and  step 
function.  In  some  intervals  they  severely  violate  the  properties  of  the  statistical  functions  they 
approximate:  the  PDF  approximation  can  be  negative  and  the  SF  approximation  can  be  negative 
or  greater  than  1.  Nevertheless,  the  amplitude  PDF  and  SF  approximations  resulting  from  (41) 
through  (43)  are  acceptable  in  that  their  relative  errors  are  less  than  -2  for  r  <  4.8  [5]. 

2.5  Approximate  Fluctuation  Models  Based  on  Discrete  Rayleigh  Mixtures 

A  formula  for  the  characteristic  SF,  be  it  exact  or  approximate,  can  be  used  to  find  approximate 
discrete  Rayleigh  mixture  formulas  for  the  statistical  functions  of  the  Rayleigh  parameter  and 
amplitude.  In  particular,  note  that 


^  l+fc  «,(rfc  -  r)  (44) 

and 


W,= 


rirk-\-l+rk)/2 

A'^k-\-rk~.l)/2 


fr(T)  dr  = 


Sr 


fTk  + 

V  2 


Tk+l  + 

2  J 


(45) 
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Rayleigh  Parameter  (x) 


Fig,  7  —  The  exact  and  approximate  forms  of  the  characteristic  SF  for  the  Rayleigh  clutter  fluctuation  model. 
The  exact  form  is  a  reversed  unit  step  function.  The  approximation  is  based  on  the  Fade  approximation  to  the 
amplitude  CF, 

where  To  =  0  and  the  last  weight,  W^,  is  calculated  using  =  Tw  +  —  'rjvr-i)/2.  It  follows 

that 

(46) 

ik=l 

Since  the  Rayleigh  parameter  axis  can  only  be  sampled  at  a  finite  number  of  points,  and  because  an 
approximation  to  the  characteristic  SF  may  be  used,  it  Wcis  found  that  to  ensure  the  weights  sum 
to  one,  any  negative  weight  should  be  eliminated  and  the  remaining  ones  should  be  ‘^normalized” 
via  the  mapping 


Wt 

Wk<- 

"  w,  +  ...  +  w^ 

(47) 

Substituting  (46)  into  (6)  through  (8)  yields 

(48) 

Snir)  «  1]  Wfc  exp 

(49) 

^(/>)  ~  exp  ^  )  . 

(50) 

Consider  the  fluctuation  models  presented  in  the  previous  section.  For  the  Nakagami  fluctuation 
model  in  (19)  through  (21)  with  N  —  0.5  and  =  1,  inversion  of  the  approximation  of  the  CF  in 
(24)  (as  a  Hankel  transform)  leads  to  [5] 

fR{r)  ^  2  Re{(-19.78778  + j258.7502)rAo  ((2.804626 +  j0.2884542)r) 

+  (30,66573  -  jl02,2120)  rKo  ((2.757823  +  j0.8901511)  r) 

+  (-8.790265  +  ;9.696267)  rKo  ((2.647426  +  ;1.590556)  r)}. 


(51) 
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Fig.  8  ~  The  relative  error  between  the  exact  form  of  the  Nakagami  amplitude  PDF  and  a  discrete  Rayleigh  mixture 
derived  from  a  Pade-based  approximation  of  the  characteristic  PDF  (solid  line),  a  discrete  Rayleigh  mixture  derived 
from  the  exact  form  of  the  characteristic  PDF  (dotted  line),  and  the  Pade-based  approximation  to  the  amplitude 
PDF  (dashed  line). 

The  corresponding  approximation  to  the  amplitude  SF  is 

SR{r)  ^  2  Re{(2.407845  +  j92.01074)  rK^  ((2.804626  +  jO.2884542)  r) 

+  (-0.76368991  -  j36.81605)  rK^  ((2.757823  +  j0.8901511)  r) 

+  (-0.8228678  +  y4.156900)  tKx  ((2.647426  +  jl.590556)  r)}.  (52) 

Approximate  discrete  Rayleigh  mixtures  to  the  amplitude  PDF  and  SF  were  found  from  the  exact 
forms  of  the  characteristic  PDF  and  SF  in  (22)  and  (23),  as  well  as  from  their  Pade-based  approx¬ 
imations  in  (25)  and  (26).  Sampling  of  the  Rayleigh  parameter  axis  was  accomplished  using  the 
formula 


T-fc  =  (  -  )  k  =  l,.  (53) 

\  '^max  J 

with  =  2,  Ar  =  r^axlkmax^  ^max  “  30  and  6  =  3.  This  produces  a  sampling  of  the  Rayleigh 
parameter  axis  that  is  ^^finer”  near  r  =  0,  This  was  done  to  derive  accurate  discrete  Rayleigh 
mixture  approximations  to  the  amplitude  PDF  and  SF  near  r  =  0.  For  this  fluctuation  model, 
linear  sampling  (c  =  1)  produces  less  accurate  approximations. 

Figures  8  and  9  show  the  relative  errors  of  the  amplitude  PDF  and  SF  approximations  when 
compared  to  the  closed  forms  in  (19)  and  (20).  Both  types  of  approximation  are  accurate  for  r  <  2.5 
where  the  relative  error  is  generally  less  than  -2. 

For  the  K-type  fluctuation  model  in  (28)  through  (32)  with  u  =  0.5  and  b  =  \/3,  inversion  of 
the  approximation  of  the  CF  in  (33)  (as  a  Hankel  transform)  leads  to  [5] 


/MO  ^  10.04852  r Ao(F712555r)  -  7.166959  rAo(1.844955r) 

-1.600556  rA"o(2.411114r)  -  0.9000533  rAo(4.455323r) . 


(54) 


Target  Detection  Using  a  Linear  Sum  of  Matched  Filter  Outputs 


13 


Fig.  9  —  The  relative  error  between  the  exact  form  of  the  Nakagami  amplitude  SF  and  a  discrete  Rayleigh  mixture 
derived  from  a  Pade-based  approximation  of  the  characteristic  SF  (solid  line),  a  discrete  Rayleigh  mixture  derived 
from  the  exact  form  of  the  characteristic  SF  (dotted  line),  and  the  Pade-based  approximation  to  the  amplitude 
SF  (dashed  line). 

The  corresponding  approximation  to  the  amplitude  SF  is 

SR{r)  5,867560  riri(L712555r)  --  3.884627  riri(L844955r) 

-0.6638242  riri(2.411114r)  -  0.2020175  riri(4,455323r) .  (55) 

Approximate  discrete  Rayleigh  mixtures  to  the  amplitude  PDF  and  SF  were  found  from  the  exact 
forms  of  the  characteristic  PDF  and  SF  in  (31)  and  (32),  as  well  as  from  their  Pade-bsised  approx¬ 
imations  in  (34)  and  (35).  Sampling  of  the  Rayleigh  parameter  axis  was  accomplished  using  (53) 
with  Trnax  =  6,  Ar  =  r^axl^max,  ^max  =  30  aud  £  =  2.  Figures  10  and  11  show  the  relative  errors  of 
the  PDF  and  SF  approximations  when  compared  to  the  closed  forms  in  (28)  and  (29). 

Generally  the  ^tails”  of  the  discrete  Rayleigh  mixture  approximations  to  the  PDF  and  SF  are 
not  as  accurate  as  those  derived  from  the  Pade  approximation  to  the  CF.  With  proper  sampling  of 
the  Rayleigh  parameter  axis,  the  peaks  of  the  Rayleigh  mixture  approximations  to  the  amplitude 
SF  are  accurate  approximations  to  the  true  SF  and  thus  provide  accurate  approximations  to  useful 
target  detection  probabilities,  specifically,  0.1  <  <  0.9. 

3.  MODELING  MULTIPLE  OBSERVATIONS 

3.1  The  Observation  Vector 

Multiple  observations  of  a  target  can  be  acquired  by  processing  echoes  generated  from  time- 
sequential  transmissions,  frequency  sub-bands  of  a  single  echo,  or  a  single  echo  steered  to  different 
bearing  angles.  Suppose  that  Nd  complex  matched  filter  outputs  are  available  and  are  denoted  by 

2/fc  =  +  JVk  for  k  =  .  (56) 

These  complex  numbers  can  be  arranged  into  a  vector  of  length  N^: 
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Fig.  10  -  The  relative  error  between  the  exact  form  of  the  K-type  amplitude  PDF  and  a  discrete  Rayleigh  mixture 
derived  from  a  Pade-based  approximation  of  the  characteristic  PDF  (solid  line),  a  discrete  Rayleigh  mixture  derived 
from  the  exact  form  of  the  characteristic  PDF  (dotted  line),  and  the  Pade-based  approximation  to  the  amplitude 
PDF  (dashed  line). 


Fig.  11  -  The  relative  error  between  the  exact  form  of  the  K-type  amplitude  SF  and  a  discrete  Rayleigh  mixture 
derived  from  a  Pade-based  approximation  of  the  characteristic  SF  (solid  line),  a  discrete  Rayleigh  mixture  derived 
from  the  exact  form  of  the  characteristic  SF  (dotted  line),  and  the  Pade-based  approximation  to  the  amplitude 
SF  (dashed  line). 
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y  =  [ui+JVi  ...  Unj+  ]  . 

Alternately,  they  may  be  arranged  into  a  real  vector  of  length  2iV<j; 

X  =  [  «i  Vi  ...  ]  • 


(57) 


(58) 


3.2  The  Covariance  Matrix 

In  Section  2.1  the  assumption  of  circular  symmetry  was  imposed  because  of  the  complete 
uncertainty  in  the  phase  of  the  complex  matched  filter  output.  In  the  case  of  multiple  observations, 
the  circular  symmetry  is  also  imposed  for  each  complex  matched  filter  output,  which  implies  that  the 
real  and  imaginary  parts  are  zero-mean  random  variables.  This  alone  is  not  enough  to  completely 
specify  the  joint  PDF,  for  it  is  also  necessary  to  assume  a  certain  form  of  the  cross-covariance 
between  the  real  and  imaginary  parts  of  the  matched  filter  outputs. 

In  an  echolocation  system,  it  is  reasonable  to  assume  that  the  correlations  between  the  matched 
filter  outputs  do  not  change  if  all  complex  matched  filter  outputs  experience  the  same  phase  shift. 
Consequently,  it  is  shown  in  Appendix  A  that  for  an  echolocation  system  the  2A ^  X  2,N j  covariance 
matrix  for  the  real  vector  x  has  the  special  form 

i;{xx"’}  =  R-[R,*]  (59) 

where  the  right  side  of  (59)  denotes  the  matrix  as  an  array  of  2  X  2  sub-matrices  of  the  form 


where  and  are  real  numbers.  Note  that  for  i  =  k  the  sub-matrix  is  diagonal,  which  expresses 
the  circular  symmetry  of  the  individual  complex  matched  filter  outputs.  For  t  5^  fc,  the  matrix 
expresses  the  additional  assumption  necessary  for  echolocation  problem:  antisymmetric  correlation 
between  the  real  and  imaginary  parts. 

From  the  sub-matrix  in  (60)  it  follows  that 

IPS  1  _  /  {or  i  =  k 

{y.yki  -  I  2<7i(7k{aik  -1-  J0ik)  for  i  7^  k 

which  are  the  elements  of  the  X  covariance  matrix 

^{yy"}-P 

where  denotes  the  Hermitian  transpose  of  the  complex  vector  y. 


(61) 


(62) 


Rtfc  — 


E{uiUk} 

E{viUk} 


E{uiVk} 

E  {u<Ufc} 


<Tj<7fc 


1  0 
0  1 

Pik 


~Pik 


for  i  =  k 
for  t  7^  A: 


(60) 
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3,3  The  Multidimensional  Gaussian  Model 


It  is  well  known  that  the  multidimensional  Gaussian  PDF  is  given  by 

exp 

^  (27r)^dVdit(5f  ■ 

The  CF  associated  with  this  multidimensional  PDF  is  its  2iVrf-dimensional  Fourier  transform  with 
respect  to  the  vector  of  frequency  variables 


M  —  [  ^1  ] 


(64) 

where  is  associated  with  u*  (real  part)  and  is  associated  with  (imaginary  part).  It  is  well 
known  that  the  CF  is  given  by 


$(/z)  =  exp  . 


(65) 


If  the  matched  filter  outputs  axe  arranged  as  a  complex  random  vector  whose  covariance  matrix 
obeys  the  properties  described  in  the  previous  section,  then  [9,19] 


^  exp(-y"P-^y) 

•'YtO'y  ^>^rfdet(P)  ■ 

Consider  the  random  variable  q  =  y^y.  Its  moment  generating  function  is  given  by 

E  {exp  (ty"y)}  =  E  {exp  (t^)} 

It  is  straightforward  to  show  that  [17,19] 

{exp  (ly'y)}  =  =  81^  =  ! 


(66) 

(67) 

(68) 


where  Ai  —  are  the  eigenvalues  of  the  covariance  matrix  P  and  the  last  equality  is  the  result 
of  a  partial  fraction  expansion  with  parameters  D.-  and  n,-.  Inverting  this  function  (as  a  Laplace 
transform)  yields  the  PDF  of  ^  to  be 


(69) 


The  amplitude  of  the  vector  y  is  w  =  so  by  a  transformation  of  random  variable,  it  follows 
that 


...  ^  2Di  2n-l  ( 

'  “’’vr)' 

Integration  of  (70)  yields 

If  Ui  —  1  for  all  i  (which  is  typically  the  case),  then  (71)  reduces  to 

Sw{w)  ~  ^  A  exp  (^-^)  • 


(70) 


(71) 


(72) 
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3.4  The  Spherically  Invariant  Random  Vector  Model 


A  generalization  of  a  random  process  that  follows  the  Rayleigh  mixture  fluctuation  model  in 
Section  2.2  is  the  spherically  invariant  random  process  (SIRP)  that  produces  spherically  invariant 
random  vectors  (SIRVs)  [1,2,3,8,13-15,19,20].  The  PDFs  for  SIRVs  are  a  generalization  of  (6)  and 
are  given  by 


/x(x)  = 


1  1 


(73) 


The  CF  for  a  SIRV  is  a  generalization  of  (8)  given  by 

$(/i)  =  $  exp  /r(r)  dr. 


(74) 


(75) 


It  can  be  show  that  the  PDF  for  a  SIRV  has  an  alternate  representation  that  is  a  generalization 
of  the  Hankel  transform.  In  Reference  13  it  is  shown  that 


h(x^R-»x) 

’  (27r)^Vdet(R) 


(76) 


/Y(y)  = 


h  (2y"P“^y) 


z-^<<det(P) 

where  h(*)  characterizes  the  square  magnitude  of  the  vector  y  and  is  given  by 


(77) 


'*  (C)  =  (v^/>)  dp-  (78) 

Note  that  (77)  is  a  generalization  of  the  multidimensional  Gaussian  PDF  and  reduces  to  this  PDF 
when  h(Q  =  exp(— C/2).  Furthermore,  (76)  and  (77)  show  that  .^ynp-iy  ^  ^x^R“^x/2  [9]. 

3,5  Approximate  Amplitude  Fluctuation  Models  Based  on  Fade  Approximations 


The  approximations  to  the  PDF,  SF  and  CF  for  single  observations  that  are  based  on  Pade 
approximations  to  the  CF  (presented  in  Section  2.3)  can  be  immediately  generalized  to  the  case 
of  multiple  observations.  From  (78)  it  can  be  seen  that  h(-)  is  only  dependent  upon  the  circularly 
symmetric  CF,  which  is  given  by  (75)  and  is  approximated  by 


M 


_ ^ _ 


(79) 


where  division  of  the  argument  by  2rj  maintains  the  expected  power  of  each  element  of  x  as  given 
by  the  matrix  R.  It  should  be  noted  that  77  is  not  a  free  parameter  but  is  given  by 


-/=E 


(80) 
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and  constrains  the  second  radial  moment  to  be  E  {x^R“‘x}  =  2Nd.  It  is  shown  in  Appendix  B 
that  T]  is  proportional  to  the  expected  power  of  the  matched  filter  output  for  a  single  observation. 
Formal  inversion  of  (79)  (as  a  Hankel  transform)  using  Eq.  (2)  on  page  434  of  Reference  18  yields 


M 

k=l 


2’"*-‘r(mfc) 


W-^vPkQ 


It  follows  from  (77)  that 


/Y(y) « 


(27])^“  f.  A,  iV-4my«P-^y) 

jr"ddet(P);^  2’”*:-‘r(mfc)  (v'-dT/p^y^P-iy)^''"’"* 


(81) 


(82) 


An  alternate  way  to  arrive  at  (82)  is  to  substitute  the  Pade-based  approximation  to  the  character¬ 
istic  PDF  in  (16)  into  (74)  and  carry  out  the  integration  using  Eq.  (15)  on  page  183  of  Reference 
18.  This  avoids  an  apparent  problem  with  the  Hankel  transform  inversion  of  (79):  the  inversion 
integral  as  given  by  Eq.  (2)  on  page  434  of  Reference  18  is  valid  only  when  0  <  Nd<  2mk  —  3/2,  a 
restriction  that  is  easily  violated  in  practice"*. 


Invoking  the  multidimensional  transformation  of  random  variable  r  =  y'y^p-iy  yields  [13] 


Ur) 


^  y 


r(TOt) 


(\/-47?Pfcr)  . 


(83) 


Closed-form  integration  of  (83)  is  only  possible  when  m*  =  1  (which  is  usually  the  case)  or  =  1, 
thus 


~  ^  Aki-pu)'"'^'^  (\/-4mr)  -  (84) 

For  m*  >  1,  approximations  to  the  integral  of  (\/— 47ip*r)  must  be  used,  some  of 

which  may  be  found  in  Reference  11.  For  example,  by  using  the  large- argument  asymptotic  form 

^^.(0«V^exp(-t)  (85) 

in  (83)  and  integrating,  it  follows  that 


It  is  demonstrated  in  Reference  5  that  the  Pade-based  approximations  are  generally  suitable  for 
approximating  the  tails  of  the  PDF  and  SF  in  single-observation  detection  problems.  This  allows 
accurate  calculation  of  false  alarm  probabilities,  which  are  typically  less  than  0.001.  It  is  expected 
that  this  is  also  true  for  multiple-observation  detection  problems. 

^We  are  reporting  the  work  of  earlier  researchers  in  this  field  who  ignored  this  restriction. 
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3.6  Approximate  Fluctuation  Models  Based  on  Discrete  Rayleigh  Mixtures 


Generalization  of  the  discrete  Rayleigh  mixture  statistical  functions  in  Section  2.5  is  also  im¬ 
mediate.  Thus,  it  follows  that  if  the  characteristic  PDF  is  given  by  (46),  then  discrete  Rayleigh 
mixture  approximations  to  (73),  (74)  and  (75)  are 


N 


/x(x) 


exp 


2 


(TfcR)' 


A/det(rfcR) 


(87) 


/v(y)  ~  ^ 


(88) 


1=1 


(89) 


Since  the  individual  terms  composing  (87)  and  (88)  are  Gaussian  PDFs,  the  PDF  and  SF  of  the 
amplitude  w  =  \/y^y  =  \/x^x  are  simply  mixtures  based  on  (70)  and  (71).  Thus, 


2D. 


Jb=l 


iS  (T.-'0-‘r(n,) 


2n-l  f 

exp  ( - ^ 

\ 


(90) 


k=i 


(91) 


where  A,,  n,  and  Di  are  as  defined  in  (68).  If  n,-  =  1  for  all  i  (which  is  typically  the  case),  then  (91) 
reduces  to 


it=i 


(92) 


3.7  The  Receiver  Structure 

The  receiver  structure  for  multiple  observations  is  shown  in  Figure  12  and  is  a  generalization 
of  the  normalized  matched  filter  employed  for  processing  a  single  echo.  The  normalization  is 
accomplished  by  the  multiplication  of  the  observation  vector,  y,  by  the  inverse  square  root  of  the 
clutter  covariance  matrix,  Pc,  where  Pc  =  Pc  ^Pc  The  resulting  normalized  observation  vector, 
z,  whose  covariance  matrix  is  the  identity  matrix  when  the  input  is  only  clutter,  is  followed  by  an 
envelope  detector.  If  the  detection  statistic  ,r,  exceeds  the  threshold,  7,  a  target  is  declared  present 
(hypothesis  Hi)^  otherwise,  only  clutter  is  assumed  to  be  present  (hypothesis  Hq). 

If  the  clutter  and  target  fluctuation  models  follow  any  of  the  classic  Swerling  models,  then 
the  receiver  structure  shown  in  Figure  1  is  ^^optimal”  in  the  sense  that  it  maximizes  a  likelihood 
function  [4].  The  receiver  structure  is  not  optimal  in  any  sense  for  detecting  targets  in  non- 
Gaussian  clutter.  It  is  addressed  here  because  it  is  a  simple,  well  understood  structure  that  is  easy 
to  implement.  Furthermore,  unlike  receivers  based  on  the  likelihood  ratio,  it  is  not  dependent  upon 
the  statistical  nature  of  the  target  or  clutter. 
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Fig.  12  -  The  multi-dimensional  normalized  matched  filter  receiver  structure. 

3.8  Setting  the  Detection  Threshold 

The  detection  threshold,  7,  is  set  by  solving  the  non-linear  equation 


5'«(7)-P/.  =  0 


(93) 


where  Sji(r)  is  either  an  exact  or  approximate  form.  Using  Newton’s  method  results  in  the  following 
recursive  equation  that  usually  converges  to  the  solution  of  (93): 


=l(k  + 


Pfa  — 

fnilk) 


(94) 


where  70  is  an  initial  value  (guess)  of  the  threshold. 


3,9  Calculating  the  Probability  of  Detection 


The  Rayleigh  mixture  approximations  can  be  used  to  calculate  detection  probabilities.  Consider 
the  following  approximations  to  the  CFs  for  the  target  and  clutter  amplitude  fluctuation  models: 

Nt 

$(/x|target)  fa  exp  ( — 

$(^|clutter)  fa  ^^Wc,i  exp  .  (96) 


Since  the  CFs  are  multidimensional  Fourier  transforms  of  the  PDFs,  the  CF  for  the  target-plus- 
clutter  case  can  be  found  by  multiplying  the  CFs: 


=  $(/^|target)$(/x|clutter) 

Nt,Nc 

jt=i,/=i 


(■^r.fcR-r  +  'rc.iR'c)^* 


(97) 


Since  the  term  in  the  summation  above  are  merely  the  CFs  associated  with  multidimensional 
Gaussian  PDFs,  from  (63)  through  (66)  it  follows  that  the  PDFs  associated  with  this  CF  are 
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/x(x|^t) 


(27r)'''<< 


exp  f-|x^  (rr.fcRr  +  Tc.jRc)  '  x] 
5-  E  Wr.,Wcj-  '■  - 


■^det(r7’ fcRx  +  rc,fRc) 


(98) 


1  exphy^(r^,,P^  +  rc,Pc)-^y] 


Nt.Nc 

/v(y|^l)^4:  E  Wr,.Wc, 


k=l,l^l 


det(rx^fcPx  +  Tc^tPc) 


(99) 


Detection  is  performed  on  the  magnitude  of  the  normalized  multi-dimensional  output,  z  —  Pc^^^y, 
Thus,  by  this  transformation  of  random  variable,  the  PDF  of  z  is 


1 

A(z|Ra)«—  5:  Wr.,Wc, 


exp 

-z" 

det  1 

1 

(100) 


The  PDF  of  r  =  y^y^p-iy  =  y/z^z  =  tv  follows  from  (70)  and  is 


2Di,jk,z 


(101) 


[tr  A”lt'r  («,,*,) 

where  Xi^k,i  *  -  -  Si-re  the  eigenvalues  of  the  covariance  matrix  Tc,zI  +  Tr^^Pc^^^P^Pc^^^  and 

ni^k,i  •  •  • » denote  their  respective  powers  in  the  partial  fraction  expansion.  Integration  of 
(101)  leads  to 


Nj>,Nc 

SRirlHi)^  Wr^kWc, 

fc=i,i=i 


^  A,M 

“r(n,.fc,j)  V  *’*■*’ A, 


If  =  1  (which  is  usually  the  case),  then  (102)  simplified  to 

'^d  f  r'^  \ 


N^,Nc 

SnirlHi)^  E 

k=lA=l 


(102) 


(103) 


The  SFs  in  (102)  and  (103)  provide  the  detection  probabilities  when  evaluated  at  the  detection 
threshold  (r  =  7). 


3.10  Examples 

The  Accuracy  of  Target-Plus-Clutter  Models 

The  example  in  this  section  quantifies  the  possible  accuracy  of  using  Rayleigh  mixture  ap¬ 
proximations  to  calculate  the  statistical  functions  of  a  non-Rayleigh  target  in  non- Rayleigh  clutter 
for  multiple  observations.  A  case  is  examined  where  the  individual  multidimensional  statistical 
functions  are  known  in  closed  form  but  the  PDF  for  the  sum  of  random  vectors  is  not.  In  general 
the  statistical  functions  of  the  target-plus-clutter  case  cannot  be  found  in  closed  form.  However, 
it  is  possible  to  compare  the  one-dimensional  (complex)  projection  of  the  joint  PDF  with  their 
corresponding  exact  analytic  forms  in  special  cases. 

Consider  the  case  of  four  complex  matched  filter  outputs  (N^  =  4)  for  a  K-type  target^  (yi)  in 
K-type  clutter  (y2),  each  with  joint  PDFs  between  their  real  and  imaginary  parts  given  by 

^Because  the  resulting  target-plus-clutter  amplitude  PDF  and  SF  can  be  expressed  in  closed  form  (shown  later), 
this  choice  of  target  model  is  only  a  mathematical  convenience. 
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=  3(4,).1(M,')*'"  (b.^y^T'y)  fori  =1.2 


where  —  0.5,  62  =  0.75  and 


M,= 


0.7092 

0.5390  -  0.0385; 
0.4749  +  0.0186; 
0.5392  -  0.0669; 


0.5390  +  0.0385; 
0.6577 

0.4587  +  0.0437; 
0.5092  -  0.0090; 


0.4749-0.0186; 
0.4587  -  0.0437; 
0.6170 

0.4891-0.0391; 


0.5392+0.0669; 
0.5092  +  0.0090; 
0.4891  +  0.0391; 
0.6930 


M2  = 


0.6532 

0.4995  +  0.0519; 
0.4764-0.0127; 
0.4971  +  0.0426; 


0.4995  -  0.0519; 
0.6559 

0.4902  -  0.0463; 
0.5167  +  0.0137; 


0.4764  +  0.0127; 
0.4902  +  0.0463; 
0.6168 

0.4815  +  0.0378; 


0.4971-0.0426; 
0.5167-0.0137; 
0.4815  -  0.0378; 
0.6701 


The  characteristic  functions  of  these  models  are 


$(p)  = 


1 

l  +  ip/bif  ’ 


for  i  =  1, 2. 


(104) 


(105) 


(106) 


(107) 


The  CF  for  the  new  random  variable  y  =  yi  +  y2  is  the  product  of  the  CFs  in  (107)  whose  partial 
fraction  expansion  is 


m  = 


blbl  r  1 

bi-bi  [i  +  ip/biY 


1 

1 + (p/b2y. 


(108) 


which  can  be  inverted  (as  a  Hankel  transform)  to  yield  the  one-dimensional  amplitude  PDF  and 
SF  as 


/s(^) 

5'ii(r) 


blbl 

bl-bl 

b\bl 


r[Ko{b,r)-Ko(b,r)] 


bl-bf 


(109) 

(110) 


As  noted  above,  the  closed-form  expressions  for  the  PDF  and  SF  of  x  cannot  be  found,  but 
the  projection  of  the  random  variables  given  by 

r  =  Vx^Ax  (111) 

where 

A  =  Diag[l,0,0,0]  (112) 

or  any  other  permutation  of  the  elements  of  A,  foUow  the  statistical  functions  in  (109)  and  (110). 

Rayleigh  mixture  approximations  to  (109)  and  (110)  were  calculated  using  the  procedure  in 
Section  3.9  using  the  sampling  of  the  Rayleigh  parameter  axis  given  by  (53)  with  =  15/^^^, 
kmax  —  50,  At  =  r^ax/fcmax?  ^^d  6  =  2.  Thc  relative  errors  of  the  Rayleigh  mixtures  are  shown  in 
Figure  13.  Except  for  a  small  interval  near  r  =  0,  these  approximations  show  good  agreement  with 
the  closed  form  expressions. 
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Fig.  13  -  The  relative  errors  between  the  closed  form  SF  and  PDF  of  the  projection  and 
their  Rayleigh  mixture  approximations. 


Calculations  of  ROC  curves 


Consider  the  case  of  a  Nakagami  fluctuating  target  in  clutter  following  a  K-type  fluctuation 
model.  The  parameters  of  these  models  are  the  same  as  those  presented  in  Section  2.4.  Let  Nd  =  4 
with  the  elements  of  the  complex  covariance  matrix  Px  following  the  Markov  form 

=  (113) 

Let  the  complex  covariance  matrix  Pc  be  the  identity  matrix,  thus  modeling  the  case  of  uncorrelated 
clutter.  Finally,  let  the  SCR  be  defined  by 


trace  (Pc) 


(114) 


which  is  the  total  target  echo  power  (from  all  matched  filter  outputs)  divided  by  the  total  clutter 
power. 


The  receiver  operating  characteristic  (ROC)  curves  were  calculated  for  this  combination  of 
fluctuation  models  using  the  procedure  outlined  Sections  3.8  and  3.9.  Pade-based  approximations 
to  the  SF  (described  in  Section  3.5)  were  used  to  calculate  the  detection  threshold.  Discrete 
Rayleigh  mixture  models  (described  in  Section  3.6)  were  used  to  calculate  P^.  In  particular,  the 
exact  form  of  the  characteristic  SF  for  the  Nakagami  model  was  used  to  calculate  its  approximate 
discrete  Rayleigh  mixture  model,  whereas  the  Pade-based  SF  approximation  was  used  for  the  K- 
type  approximate  mixture  model. 


The  ROC  curves  plotted  on  total  probability  scales  are  shown  in  Figures  14  through  16,  each 
showing  Prf  for  a  different  Four  ROC  curves  are  shown  in  each  figure:  Rayleigh  target  in 

Rayleigh  clutter,  =  1  (dotted  line);  Nakagami  target  in  K-type  clutter,  N4  —  4,  uncorrelated 
target  (dashed  line);  Nakagami  target  in  K-type  clutter,  =  4,  correlated  target  (solid  line); 
Nakagami  target  in  K-type  clutter,  Na  =  I  (dot-dash  line). 
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Fig.  14  -  The  ROC  curves  for  P/*  =  0,001:  Rayleigh  target  in  Rayleigh  clutter,  Nd  =  1  (dotted  line);  Nakagami 
target  in  K-type  clutter,  Nd  =  4,  uncorrelated  target  (dashed  line);  Nakagami  target  in  K-type  clutter,  Nd  =  4, 
correlated  target  (solid  line);  Nakagami  target  in  K-type  clutter,  Nd  ^  I  (dot-dash  line). 

As  expected,  using  multiple  observations  yields  higher  detection  probabilities.  The  curves 
for  the  uncorrelated  target  are  nearly  identical  in  shape  to  the  ROC  curve  for  correlated  target 
observations,  but  are  shifted  approximately  1  dB  to  the  left.  The  ROC  curves  for  the  Rayleigh 
target  in  Rayleigh  clutter  are  well  known  and  obey  the  equation 

P<1  =  (115) 

They  have  been  included  for  reference  since  they  may  be  regarded  as  standard  set  of  ROC  curves. 
Clearly  the  mathematical  form  of  the  fluctuation  models  have  a  profound  effect  on  the  performance 
of  a  detection  system. 

It  was  found  that  when  using  total  probability  scales  to  plot  ROC  curves,  as  was  done  in  Figures 
14  through  16,  small  variations  in  the  curves  arising  from  the  use  of  different  approximations  will 
be  exaggerated  for  >  0.9,  This  should  not  be  of  great  concern  because  it  is  the  portion  of  a 
ROC  curve  where  0.1  <  P^  <  0.9  that  is  the  most  useful. 

4.  CONCLUSIONS 

The  mathematical  techniques  presented  in  this  report  allow  accurate  computation  of  detection 
probabilities  over  a  range  of  values  considered  critical  in  the  performance  evaluation  of  echoloca- 
tion  systems.  The  Pade-based  approximations,  which  generally  provide  good  approximations  to 
PDF  tails,  are  suitable  for  calculating  detection  thresholds.  The  approximations  based  on  discrete 
Rayleigh  mixtures  provide  good  approximations  to  PDF  modes  (peaks)  that  generally  affect  the 
probability  of  detection. 

Though  the  techniques  presented  here  appear  complicated,  especially  in  light  of  the  matrix 
calculations  that  must  be  performed,  it  was  found  that  implementation  was  straightforward  using 
Matlab. 
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Fig.  15  -  The  ROC  curves  for  JV.  =  0.0001:  Rayleigh  target  in  Rayleigh  clutter,  ==  1  (dotted  line);  Nakagami 
target  in  K-type  clutter,  Nd  =  4,  uncorrelated  target  (dashed  line);  Nakagami  target  in  K-type  clutter,  JV<|  =  4, 
correlated  target  (solid  line);  Nakagami  target  in  K«type  clutter,  Nd^X  (dot-dash  Une). 


Fig.  16  -  The  ROC  curves  for  P/a  =  0.00001:  Rayleigh  target  in  Rayleigh  clutter,  Nd  =^X  (dotted  line);  Nakagami 
target  in  K-type  clutter,  Nd  =  4,  uncorrelated  target  (dashed  line);  Nakagami  target  in  K-type  clutter,  Nd  =  4, 
correlated  target  (solid  line);  Nakagami  target  in  K-type  clutter,  Nd  =-X  (dot-dash  line). 
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Implementation  of  the  techniques  presented  here  requires  some  human  intervention.  For  ex¬ 
ample,  choosing  the  order  of  the  Fade  approximations,  the  sampling  of  the  characteristic  SF  and 
assessing  the  quality  of  the  resulting  statistical  function  approximations  requires  some  judgment  by 
an  analyst.  Consequently,  the  techniques  presented  herein  are  best  suited  for  off-line  performance 
prediction. 
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Appendix  A 

THE  COVARIANCE  MATRIX  FOR  ECHOLOCATION  PROBLEMS 


When  two  or  more  complex  matched  filter  outputs  are  used  to  detect  a  target,  it  is  possible 
that  an  identical  phase  shift  can  appear  across  all  the  of  the  outputs.  That  is,  if  the  fc-th  matched 
filter  output  ttfc  +  jVk  then 

«fc  +  JVk  +  jv*)  expOV’)  =  [«*  cos(V>)  -  Vk  sin(^)]  +  jK  sin(V’)  +  cos(^)].  (Al) 

This  can  occur  in  the  following  situations: 

•  If  observations  are  obtained  from  a  single  matched  filter  output  that  is  sampled  in  range  delay 
or  is  sampled  at  the  same  range  delay  over  several  pings,  a  common  phase  shift  can  appear  in 
all  the  observations  if  the  target  exhibits  a  slight  change  of  range  (smaller  than  a  resolution 
ceU)  between  scans^. 

•  In  frequency  diversity  system  where  observations  are  drawn  simxiltaneously  form  several  differ¬ 
ent  frequency  bands,  slight  changes  in  target  range  (again,  between  scans)  will  induce  a  phase 
change  that  wiU  be  nearly  the  same  across  aU  frequency  bands  if  the  percent  difference  between 
the  largest  and  smallest  band  carrier  frequency  is  small. 

•  If  observations  are  drawn  from  a  compact  set  of  cells  in  the  range-azimuth  display,  then  a  small, 
common  shift  in  the  steered  beam  directions  (through  reprocessing  of  the  echoes)  wiU  induce 
an  identical  phase  shift  in  the  values  of  the  sampled  matched  filter  outputs. 

In  all  of  these  situations,  it  is  reasonable  to  assume  that  the  statistical  nature  of  the  matched  filter 
outputs  do  not  significantly  change,  the  result  being  that  the  correlation  matrices,  P  and  R  as 
defined  in  Section  3.2,  do  not  change.  This  is  a  generalization  of  the  assumed  circular  symmetry 
of  the  joint  PDF  of  the  real  and  imaginary  parts  of  the  individual  complex  matched  filter  output. 

From  (Al),  the  invariance  of  the  correlation  matrix  under  the  abovementioned  phase  shift 
implies  that 

E{uiUk}  E{uiVk} 

E{ViUk}  E{ViVk} 

where 

A  =  E  {uiUk}  cos^ii))  +  E  {v,Vk}  sin^(^)  -  (E  +  E  {u.u  J)  cos(^)  sin(^)  (A3) 

B  =  E  {u,Vk}  cos^ii))  -  E  {viUk}  sin^irp)  +  (E  -  E  cos{ip)  sm(^)  (A4) 

C  =  E  {uiUfc}  cos^(^)  -  E  +  (E  -  E  cos(^)  sin(0)  (A5) 

D  =  E  {viVk}  cos^(V^)  -I-  E  {uiUk}  sin''(^)  +  (E  {u.u J  +  E  J)  cos(V^)  sin(^).  (A6) 

®Iii  radar  parlance,  a  “scan”  refers  to  a  collection  of  observations  taken  within  a  period  of  time  when  the  target 
is  within  the  antenna  beam.  Within  a  scan  the  target  may  be  regarded  as  being  stationary  in  range. 


A  B 
C  D 


(A2) 
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These  equations  can  be  expressed  in  the  following  matrix  form: 


cos^(^)  —  1  —  cos(^)  sin(^)  —  cos(^)sin(V’)  sin^(V’) 

■  A  ■ 

■  0  ■ 

cos(^)sin(^)  cos^(^)  —  1  —  sin^(V’)  —  cos(^)  sin(^) 

B 

0 

cos(^)  sin(V’)  —  sin^(V’)  cos^(^)  —  1  —  cos(^)  sin(^) 

C 

0 

sin^(^)  cos(V’)  sin(^)  cos(^)  sin(0)  cos^(V’)  —  1 

D 

0 

This  system  of  equations  can  be  simplified  by  adding  the  first  row  to  the  last,  subtracting  the 
second  row  from  the  third,  and  multiplying  the  second  row  by  cos(V’)/  sin(V’)  and  adding  the  result 
to  the  first  row.  The  resulting  reduced  set  of  equations  are 


-10  0  1 
cos^(^)  —  cos(V’)  sin(^)  —  cos(V’)  sin(V’)  —  cos*(^) 


■  A  ' 

■  0  ' 

B 

0 

C 

0 

D 

0 

(A8) 


The  first  row  of  (A8)  implies  that  A  =  D.  Substituting  this  result  back  into  the  second  line  of  (A8) 
implies  that  B  =  —C.  Therefore, 


E{uiUk}  E{uiVk} 

'  A  -C' 

E{ViUk}  E{ViV^} 

C  A 

(A9) 


This  proves  the  form  of  the  covariance  sub-matrices  defined  in  (60). 


Appendix  B 


THE  EXPECTED  POWER  FROM  PADE  APPROXIMATIONS  TO  THE 

CHARACTERISTIC  FUNCTION 


In  this  appendix  it  is  shown  that  the  expected  power  of  a  fluctuation  model  for  N4  =  1  (single 
observation)  can  be  cadculated  from  the  coefficients  and  poles  of  the  Fade  approximation  to  the  CF 
that  is  given  by  (Sections  2.1  and  2.3) 

Hp)  =  1 - ~  ^ 


where  r  is  the  amplitude  of  the  matched  filter  output  defined  in  (1).  It  is  easy  to  show  that 


^kP"^  ,  (  4n 


{p^-VkT*"  (“pO"**  (-pO”**"^' 

Substituting  (B2)  into  (Bl)  and  grouping  the  coefficients  for  each  power  of  p  reveals  that 


(B2) 


where  t)  is  the  normalization  coefficient  defined  in  Section  3.5.  This  formula  is  always  exact  if 
moment  matching  is  used  to  find  the  coefficients  of  the  Fade  approximation  in  (Bl).  This  is  because 
the  second  moment  of  the  amplitude  PDF  is  always  used  in  the  matrix  equations  describing  the 
Fade  coefficients  ((10)  through  (12)  in  Section  2.3). 
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